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ABSTRACT 

The majority of clusters in the Universe have masses well below 10^ Mq. Hence their inte- 
grated fluxes and colors can be affected by the presence or absence of a few bright stars introduced 
by stochastic sampling of the stellar mass function. Specific methods are being developed to ex- 
tend the analysis of cluster energy distributions into the low-mass regime. In this paper, we 
apply such a method to real observations of star clusters, in the nearby spiral galaxy M83. We 
reassess the ages and masses of a sample of 1242 clusters for which UBVIHa fluxes were obtained 
from observations with the WFC3 instrument on board of the Hubble Space Telescope. Synthetic 
clusters with known properties are used to characterize the limitations of the method (valid range 
and resolution in age and mass, method artifacts). The ensemble of color predictions of the 
discrete cluster models are in good agreement with the distribution of observed colors. We em- 
phasize the important role of the Ha data in the assessment of the fraction of young objects, 
particularly in breaking the age-extinction degeneracy that hampers an analysis based on UBVI 
data only. We find the mass distribution of the cluster sample to follow a power-law of index 
—2.1 ± 0.2, and the distribution of ages a power-law of index —1.0 ± 0.2 for log(M/ Mq) > 3.5 
and ages between 10*" and 10^ yr. An extension of our main method, that makes full use of the 
probability distributions of age and mass obtained for the individual clusters of the sample, is 
explored. It produces similar power-law slopes and will deserve further investigation. Although 
the properties derived for individual clusters significantly differ from those obtained with tradi- 
tional, non-stochastic models in about 30% of the objects, the first order aspect of the age and 
mass distributions are similar to those obtained previously for this M 83 sample in the range of 
overlap of the studies. We extend the power-law description to lower masses with better mass 
and age resolution and without most of the artifacts produced by the classical method. 
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1. Introduction 

Many and possibly most stars form in clus- 
ters rather than individually. Much observa- 
tional effort has therefore been invested to de- 
termine the distributions of star cluster ages 
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ble Space Telescope and especially its latest survey 
instrument WFC3, cluster samples have been de- 
tected in external galaxies down to a regime in 
which cluster luminosities overlap those of indi- 
vidual bright stars, and hence to lower masses 
(jjohnson et al. 2011). Because cluster luminosity 



and to determine the dominant mechanisms of 



cluster formation and disruption (jKroupa fc Boilv 



(and mass) distributions rise steeply towards lower 
luminosities (and masses), these deeper surveys 
have led to significantly larger cluster samples. 
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In the low-mass regime, clusters of a given 
age and mass are predicted to display a broad 
range of integrated luminosities and colors, mostly 
as a consequence of the random sampling of 
the up per part of the stella r mass function 
fiBarbaro fc Bertellil [l977t [cirard i fc Bica 199^ 

Lancon fc Mouhcinj200o[lBruzuall2002i :rC ervino fc Luridiebal would be dispersed around the real values, 
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Assigning fundamental properties to individ- 
ual clusters and studying the global properties of 
the cluster population of a galaxy are two dif- 
ferent issues. For instance, iFouesneau fc Lancon 
(I2OIOI) showed that the mass estimates deter- 
mined for low-mass clusters with continuous mod 



20061 iDeveikis et alll2008l ). For the sake of sim- 
plicity, we will refer to stellar population model 
implementations that explicitly predict luminosity 
and color distributions as stochastic models. The 
predicted distributions depend in non-trivial ways 
on the cluster mass. Stochastic models are usu- 
ally based on Monte Carlo simulations, although 
first at tempts at analytical appro aches have been 
made (jCervino fc Luridianal 120061 ). Traditionally, 
stellar population models have been implemented 
with the assumption of a fully sampled initial 
stellar mass function (IMF): each stellar mass bin 
along an isochrone contains exactly the mass frac- 
tion prescribed by the IMF, regardless of the fact 
that this may lead to unphysical non-integer num- 
bers of stars in some of the most luminous phases 
of evolution. These implementations, which we 
will refer to as continuous models, do not predict 
distributions of cluster luminosities and colors but 
rather their mean values (and sometimes their first 
moments) . 

Monte-Carlo simulations of low-mass clusters 
have shown that the range of broad-band col- 
ors predicted with this type of stochastic model 
are in rather good agreement with observed 
ones in t he Milky Way and in the Magellanic 
Clou ds I Girardi fc Bical 1993 ; Piskunov et al 
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minosity distributions are complex in the low-mass 
regime. They display multiple peaks. The most 
probable luminosities and colors differ from the 
mean color of the distribution, and conversely the 
mean properties are unlikely to occur in any actual 
cluster. Therefore ages and masses derived from 
the analysis of cluster luminosities and colors are 
expected to depend on the type of model used (i.e. 
stochastic or continuous). For different combina- 
tions ofj^road^bcuid filters, Mafz ApcUaniz (200£j) 
and IFouesneau fc LanconI (|201(]h confirmed that 
the analysis of the integrated light of synthetic 
low-mass clusters with continuous models leads 
to biases in the assignment of ages to individual 
objects. 



but that there would not be a strong global im- 
pact on the cluster mass distribution itself. This 
work also indicated that high resolution features 
in t he cluster age distribu t ion w ould be affected, 
but Fouesneau fc LangonI ( 2010l ) did not exam- 
ine whether any of these biases would be strong 
enough to modify the global, low resolution age 
distribution of a typical cluster sample. 

Th e main purpose of the article of lFouesneau fc Lancon 
(l2010h was to introduce two analysis methods for 
the interpretation of integrated fluxes of low-mass 
star clusters, both based on a large collection of 
Monte-Carlo simulations of clusters of finite num- 
bers of stars, i.e. on stochastic models. One anal- 
ysis method is a traditional minimization. The 
second one rests on the calculation of posterior 
probability distributions in age-mass-extinction 
space. The latter refers explicitly to Bayes' theo- 
rem, in which the posterior probabilities of model 
parameters given a set of data are expressed in 
terms of their prior distributions and of the prob- 
ability distributions of observational errors. Al- 
though the results from the analysis methods 
have their own model-dependencies through the 
adopted population synthesis code and the priors, 
it is clearly a conceptual improvement to move 
away from the restricted framework of continu- 
ous population synthesis models. Studies in the 
regime well below 10^ Mq, where the literature 
on the subject does not usually venture, become 
possible. This should allow us to examine cluster 
formation and disruption histories in more detail 
in the future. 



In this paper, we apply the methods of lFouesneau fc Lancon 

(fioiOi) to the WF C3 sample of star cluste rs in M 83 
first presented bv IChandar et al. ( 2010( ). Messier 
83, the "Southern Pinwheel" galaxy, is one of the 
best nearby analogs of the Milky Way. These 
data allow us to compare standard age and mass 
estimates based on continuous models with those 
based on stochastic ones for a sample of real ob- 
jects. We also investigate how the use of stochastic 
models affects the age and mass distributions of 
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the sample as a whole. 

The rest of this paper is organized as follows. 
Section 2 summarizes the data, cluster selection 
and photometry, and Section 3 summarizes and 
extends pr edictions from th e stochastic models de- 
veloped by l^uesneau&^Ta^gon QQlfl). Section 4 
compares the luminosity and color distributions of 
the colors and model predictions. Expectations 
from artificial clusters are presented in Section 5, 
and results for real clusters in M 83 are presented 
in Section 6. Section 7 discusses several key is- 
sues, including the importance of using a narrow- 
band filter in the analysis, and a direct comparison 
of results from stochastic and continuous models. 
More details relevant to this comparison are pro- 
vided in two appendices. Section 7 also present 
preliminary exploration of an alternative determi- 
nation of the age-mass distribution of the cluster 
sample of M 83, one which uses the probability dis- 
tributions of individual clusters rather than just 
their most probable values. Section 8 presents the 
main conclusions of our work. 

2. Data and Observations 

In order to take advantage of the wavelength 
coverage of WFC3, images of M83 were taken as 
part of the ERSl program 11360 (PI: O'Connell) 
through seven broad-band filters from the UV 
to the near infrared: F225W (UV), F336W (U), 
F438W (B), F555W (V), F814W (I), FllOW 
(J), F160W (H). They cover the nucleus and the 
North- Eastern part of the galaxy. The color image 
in Fig. [1] illustrates the observed region. In addi- 
tion, narrow band filter observations were taken 
for the following emission lines: [O III] (F373N), 
H/3 (F487N), O II (F502N), Ha-|-[N II] (F657N), 
[S II] (F673N). 

In this paper, we focus only on the broad band 
UBVI photometry (HST/WFC3 passbands) and 
the narrow band Ha measurements (F657N). The 
sample contains 1242 star clusters with measure- 
ments in these five band^. They were selected 
based on a "white" light image (i.e. a co-added, 
rms-normalized combination of the U, B, V and 
I filters) , in order to give roughly equal weight to 
all wavelengths. A combination of size and con- 
centration criteria measured in the V band was 



used to sepa rate individual s tars fr om clusters, as 
described in IChandar et al.l ( 2010f ). For each ob- 
ject, circular aperture photometry with an aper- 
ture radius of 3 pixels (0.0396" pix~^) produced 
integrated flux values which were then corrected 



for fo reground galactic extinction ((Schlegel et al 



Il998l . appendix B). Measurement uncertainties are 
typically spread between 0.05 and 0.25 mag for the 
B, V and I broad bands (cf. Sect. [5]). About 
one half of the objects have uncertainties larger 
than 0.25 in U and Ha filters. The average er- 
rors increase non linearly with magnitude. In- 
strumental magnitudes were converted into the 
VEGA magnitude system by applying the follow- 
ing zero-points: F336W=23.46, F438W=24.98, 
F555W=25.81, F814W=24.67, F657N=22.35. We 
as sume a distance o f 4.5 Mpc to M83 as found 
in lThim et al. (2003), which corresponds to a dis- 
tance modulus of m — M = 28.29. 

3. Population synthesis models 

We aim at studying the age-mass distribu- 
tio n of the clusters in M 83, using the method 
of Fouesneau fc LanconI ( 2010f ). The method is 
based on a large collection of Monte-Carlo simula- 
tions of individual clusters. The synthetic clusters 
are constructed with the population synthesis code 
Pegase.2n (Fouesnea u et al. in prep.), which is 
deriv ed from Pegase ( Fioc fc Rocca-Volmerang3 
1997h . As in the original (continuous) popula- 
tion synthesis code, the underlying stellar evo- 
lution tracks are t hose of the Padova group 
( Bressan et al ] Il993h . with a simple extension 
through the thermally pulsating AGB based 
on th e prescriptions of ICroenewegen fc de Jonj 
(Il993l) . The in put stellar spe c tra ar e taken from 
the library of iLeieune et al.l (|l997l) . The stel- 
lar Initial M a ss Fun ction (IMF) is taken from 
Kroupa et al.l (1993), and extends from 0.1 to 
120 Mq. Nebular emission (lines and continuum) 
is included in the calculated spectra and broad 
band fiuxes under the assumption that no ion- 
izing photon escapes. Line ratio s are computed 



as m 



Fioc fc Rocca-Volmerangd (1997). When 



^ The cluster catalog, inch photometry and positions, can be 
retrieved at: |http: //archive . stscl edu/prepds/wf c3ers| 



extinction corrections are considered, they ar e 
based on the standard law o fICardeUi et al. I (ll989l) . 
The synthetic photometry for the artificial clus- 
ters is computed using the response curves of the 
HST/WFC3 filters. A reference spectrum of Vega 
provides zero magnitude fiuxes (lBohlinll2007h . 
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Fig. 1. — Color image of the observed region of M 83 with the WFC3 instrument on-board HST. The B-band 
(F438W) image is shown in bhie, the V-band (F555W) image in green and a combination of the I-band 
(F814W) and Ha in red. The covered area represents 2.75' x 2.75' (w 3.6 x 3.6 kpc^). 
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In the context of studying M 83, the model col- 
lection is restricted to solar meta llicity, the metal- 



licity of this galaxy according to iGil de Paz et al 
()2007l) . The simulated cluster set has been 
extende d to higher masses than available in 
iFoucsnc au fc Lancon ( 2010l ). and now covers 
masses from 10"^ to 5 x 10^ Mq and ages from 
1 Myr to 20 Gyr. With 1.4 x lO'^ individual ob- 
jects, the collection is large enough to include all 
reasonably likely cluster properties. The ages of 
the synthetic clusters are drawn from a power law 
distribution with index — 1 (equal numbers of star 
clusters per logarithmic age bin), rounded to in- 
teger multiples of 10^ yr. For practical reasons 
(the need to include more massive clusters with- 
out having to recompute prohibitive numbers of 
low-mass clusters), the mass distribution in the 
new collection of models also follows a power law 
of index — 1, instead of the pre viously adopted —2 



Fouesneau fc Lancon ( 201oH ^. 



The age and mass distributions of the synthetic 
clusters, together with the values allowed for ex- 
tinction, are the main assumptions (priors) of the 
inversion. These distributions account for two ma- 
jor qualitative trends found in star-forming galax- 
ies: low mass clusters are more numerous than 
high mass clusters, and, because of a variety of ef- 
ficient disruption mechanisms, young clusters are 
more numerous than old ones. In the future, it will 
be interesting to implement an iterative inversion, 
in which the derived cluster age and mass distribu- 
tions, corrected for selection effects, are re-injectcd 
as priors until the procedure converges. 

4. Luminosity and color distributions 



Previous work by IChandar et al. I (|20ld) shows 
that the V band luminosity function of the ob- 
served clusters can be described globally as a 
power-law with index ^ — 2 down to the complete- 
ness limit located near My = —5. For the four 
other bands used here, the turn-over points in the 
luminosity distributions are located near Mu = 
-6, Mb = -5, Mi = -6 and Mh^ = -6 (in the 
absolute Vega magnitude system) . For brevity, we 
will refer to these turn-over points as "complete- 



"The number of massive clusters in tlie current collection is 
larger than strictly necessary. The choice of the power law 
index can be optimized in the future to construct smaller 
representative cluster sets. 



ness limits". Note that Chandar et al. ( 2010l) es- 
timated that this cluster catalog contains contam- 
ination at the approximately 15% level from indi- 
vidual stars in crowded regions and from chance 
superpositions of a few stars. 

Figure [2] compares the loci of the observa- 
tions and the models in three projections of color- 
magnitude space. In the left hand panels, the data 
are shown together with the complete set of mod- 
els, which will be used to assign age, mass and ex- 
tinction estimates to each individual object. The 
majority of the observed clusters lie well within 
the regions covered by the synthetic clusters. It is 
also clear that some of the observations cannot 
be reproduced with continuous population syn- 
thesis models, even when allowing for extinction. 
The clusters with the bluest (V-I) colors are those 
whose post-main sequence happens to be under- 
populated, while the clusters with the reddest (V- 
I) colors may either happen to have more luminous 
red stars than average, or be severely reddened (or 
both). 

At very young ages the color predictions are 
sensitive to the prescription adopted for the neb- 
ular emission: the nebular fluxes added to the 
stellar fluxes produces redder WFC3 V— I colors 
(to allow comparison, the red line in Fig. [2] illus- 
trates the effect of discarding the gas contribution, 
in the case of continuous models). This leads to 
the hook-like extension seen in blue at the top of 
the upper panel of Fig. [51 composed of low-mass 
clusters (a few lO'^ M©) with a number of ionis- 
ing stars in excess of the average. They do not 
have counterparts in the observed cluster sample. 
The models in that part of the diagram represent 
only 0.1% of the models with ages below 3 Myr, 
and will therefore have little effect on any of our 
results. On the other hand, the very blue colors 
observed for a subset of obviously young clusters in 
M 83 suggest that the nebular emission included in 
our models is too high for some objects; this could 
affect age-dating, mainly between 1 and 4 Myr. 

One can conclude from the two bottom-left pan- 
els of Fig. [5] that a majority of the observed clus- 
ters have masses lower than 10'' Mq. At these 
masses, the predicted optical fluxes of star clus- 
ters are spread quite widely. For instance, the 
90 % confidence intervals for the predicted fluxes 
respectively have widths of 0.3 dex at 5 Myr in 
the V band, 0.1 dex at 50 Myr in V, and 0.5 dex 
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Fig. 2. — Color-color and color-magnitude diagrams of the observed clusters (black points) and the models. 
Magnitudes are in the Vega system, in the WFC3 filters. Solid lines show the predictions obtained with the 
corresponding "continuous" population synthesis models for 10'^, 10^, 10^, and 10^ Mq. The red line in the 
top right panel shows the effect of excluding any ionized gas emission from the predicted fluxes. Colored 
dots represent the discrete clusters available in the Monte-Carlo collection. The color-code is given in the 
bottom right panel. In the right hand panels, only discrete models brighter than the observational limits in 
UBVI-f Ha are shown. Models are not reddened in any panel, but the extinction vector for Ay = 1 is shown. 
Hatched regions represent the expected colors of individual stars. 
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at 5 Myr in the Ha filter (based on our Monte 
Carlo simulations). The distributions have com- 
plex shapes. This is the mass regime where the 
use of the method of Fouesneau fc Lan(;;on ( 2010l ) 
may be expected to produce the most significant 
changes with respect to traditional approaches. A 
main purpose of this paper is to illustrate to what 
extent this does (or does not) modify results for 
individual clusters and for the sample as a whole. 

The right hand panels illustrate some of the 
complex effects of magnitude limits. The models 
shown are those brighter than the "completeness 
limits" of the observations in all five bands. Some 
clusters are observed below these limits, because 
the luminosity distributions of the observations do 
not fall off abruptly beyond their turn-over: 10% 
of the clusters are below the turn-over of the lu- 
minosity function in V or Ha, 30 to 40% in J7, B, 
and /. More complex figures are obtained when 
reddened models are included (not shown). In the 
stochastic context, a magnitude cut does not sim- 
ply reject all the clusters of a given mass that have 
reached a critical age. On the contrary the mass- 
to-light ratios in this regime have a large range 
of possible values. At a given mass and age only 
those stochastic clusters that happen to be fainter 
than the magnitude limits will be rejected, while 
those that happen to be brighter will be kept. 
From the figures shown here one should retain that 
the effects of a stochastically sampled IMF will be 
more important for the subset of young clusters 
than for the old ones, because the latter must be 
of high mass to be detected. 

The model distributions are not and should not 
be truncated in magnitude when used to assign 
fundamental properties to individual clusters. The 
presence of low luminosity clusters in the Monte- 
Carlo sample can do no harm in the analysis of a 
cluster observation. On the other hand, the empir- 
ical sample may contain a small number of clusters 
with intrinsic luminosities below the "complete- 
ness limits" , that made it into the sample only 
thanks to favourable error bars. With a truncated 
theoretical distribution, the analysis would fail to 
provide correct properties for such objects. 

In the future however, it will be worthwhile to 
construct more sophisticated versions of the trun- 
cated model distributions in the right hand pan- 
els of Fig.[2j If the original Monte Carlo collection 
can be modified to account for all observational ef- 



fects present in the M 83 sample, the resulting syn- 
thetic color-magnitude distributions can be com- 
pared directly with the observed ones. This ap- 
proach would be a direct analog of what has slowly 
become common practice in the analysis of the re- 
solved stellar color-magnitude diagrams of nearby 
galaxies. In the latter case the aim is to determine 
star formation histories. In our case current age- 
mass distributions could be estimated, without the 
detour through the individual analysis of each ob- 
served cluster. We have not yet implemented the 
tools necessary for such a study. The observational 
effects relevant to a cluster sample are more com- 
plex than those one deals within star samples (i.e. 
variable cluster surface brightness, size, effect of 
crowding, background fluctuations. .. ). Correct- 
ing for them lies beyond the scope of this paper, 
and we will rather restrict conclusions about the 
cluster population of M83 to ranges where incom- 
pleteness is not severe. 

5. Expectations from artificial data 



Fouesneau fc Lancon ( 2010^ assessed the meth- 



ods used here to estimate cluster ages and masses. 
Using standard UBVI or UBVIK photometry and 
accounting for observational errors of 5% on the 
fluxes, they showed that it was possible to de- 
termine age distributions with a resolution in age 
much narrower than in non-stochastic studies. In 
this section, we briefly describe the behaviour of 
the analysis when using photometry in the UBVI 
and Ha filters of HST/WFC3, with error bars dis- 
tributed as in the actual M 83 cluster data (Fig. [3]) . 
Most of these uncertainties are larger than 5%. 
This experiment allows us to identify potential ar- 
tifacts, and to determine what resolution in age to 
aim for in the analysis of the empirical dataset. 
Unless otherwise stated, t he second method of 
Fouesneau fc LanconI ( 2010l ) is used in the analysis 
here and in subsequent sections, i.e. we maximize 
posterior probabilities. We refer to Sect. 17. 41 for a 
discussion of a more complete exploitation of the 
posterior probability distributions. 

The sample of synthetic clusters used as a mock 
dataset contains 1000 objects and is built as fol- 
lows: (i) the number of stars in each cluster 
is drawn from a power-law distribution with in- 
dex —2, between lO'^ and 10^ (which corresponds 
roughly to masses between 500 and 5 x 10^ Mq), 
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(ii) the logarithms of the cluster ages are drawn 
from a uniform distribution, (iii) the uncertainties 
are assigned randomly to uncertainties taken from 
the list of clusters observed in M83 (see Fig.[2|), 
and (iv) extinction is added with Ay distributed 
uniformly between and 3. Note that no noise 
is added to the synthetic fluxes (the uncertainties 
assigned are used on l y in the inversion). As in 
Fouesneau Sz LanconI (|2010f ). the physical quanti- 
ties used in the inversion are fluxes, not magni- 
tudes or colors. Errors initially provided in mag- 
nitudes are reinterpreted as symmetric errors on 
the fluxes. 

The top panel of Figure |4] illustrates the input 
sample of synthetic clusters. The color scale refers 
to the continuum-subtracted Ha emission of the 
model|3, which is directly related to the tempera- 
tures and luminosities of the hot bright stars. This 
Ha emission is significant only for clusters younger 
than 8 Myr (but some young clusters show none 
as they happen to contain no ionizing star). Note 
that unlike Fig. [2l the synthetic clusters shown 
in Fig. |4] are reddened, hence they trail down the 
reddening vector along diagonal lines. There is 
a region of overlap in broad band colors between 
models with and without Ha emission (i.e. the 
dark red points which represent clusters with no 
Ha emission overlap with the points of other col- 
ors which do emit Ha ). Including F657N data 
is expected to provide critical information in this 
regime. 

The second and third panels of Fig.|4] show 
that ages are recovered without large biases and 
that errors in the derived log(j4/?;r) have a quasi- 
normal distribution with a standard deviation of 
0.14 dex (the global offset of 0.011 dex is not signif- 
icant, it can be explained by a handful of outliers). 
Considering the substructure seen in the middle 
panel, we will not attempt to interpret features 
narrower than about 0.33 dex (FWHM of a Gaus- 
sian fit) in the age distributions obtained for the 
clusters of M 83. The equivalent figures for the in- 
put and output masses show a dispersion of about 
0.1 dex. 

Large errors in age occur for a few percent of 
the synthetic clusters, and mostly at ages of one 
or a few Gyr. The origin of these is explained in 
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Fig. 3. — Photometric uncertainty distributions 
for the 1243 M 83 cluster fluxes in the five photo- 
metric bands used in this study. 



''^The continuum subtraction is performed for this illustra- 
tion only; total narrow-band fluxes are used in the analysis. 
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iFouesneau fc Lanpon ( 201(J) (their Figs. 5, 6 and 
8). Regions of high model density in color-color 
space act as attractors in our analysis. If the lo- 
cus of dereddened versions of the photometry of 
a cluster approaches such a densely populated re- 
gion, the predominant age-mass properties of the 
models located there will be considered most prob- 
able. If the cluster observed is in fact reddening- 
free, the result will be an overestimated extinction 
together with an underestimated age (the well- 
known age-extinction degeneracy is also present 
in the stochastic context). If the observed cluster 
is truly reddened, it sometimes happens that the 
analysis underestimates extinction and therefore 
overestimates age. This last situation, however, is 
unlikely to occur above 1 Gyr in the M83 sam- 
ple, because most old reddened clusters are below 
our detection limit. A small percentage of young 
reddened clusters could be affected. 

6. M 83 cluster ages and masses 

We can now go back to M83 and estimate 
age, mass and extinction for each of the clus- 
te rs in the sample in th e stochastic context. As 



Fouesneau fc LangonI (jioiO ). the assigned esti- 



mates are the age-mass-extinction triplets that 
maximize posterior probability, i.e. maximize 
the probability of the observed set of fluxes, 
given the underlying age and mass distributions 
of the model collection. The prior distribution 
in log{age) is uniform, the prior distribution in 
log(mass) is also uniform, the metallicity is solar, 
and Ay is allowed to take any value between 
and 3 in steps of 0.2. We recall that the com- 
pleteness limits of the observations are not used 
in any manner during the determination of indi- 
vidual ages and masses. A discussion of the errors 
on the estimates as derived from posterior proba- 
bility distributions in age-mass space can be found 
in Sect.m 

The age-mass-extinction distributions resulting 
from the analysis are shown in Fig. [S] Derived 
masses range from the low limit of our model cat- 
alog (10^ M©) to about 510^5 M©. Ages are dis- 
tributed between a few Myr and about 1 Gyr, 
with only a few candidates for ages older than 
1 Gyr. The relative lack of old clusters was ex- 
pected from the detection limits (Sect. 14]). The 
middle panel of Fig. [5] shows that very young clus- 
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Fig. 4. — Analysis of the colors of the mock data 
(Sect. [5]). Top: Input sample. Colors highlight 
synthetic clusters with Ha line emission. Middle 
&: bottom: Comparison between input ages and 
derived ages. The linear regression shown has a 
slope of 0.971. 
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Fig. 5. — AKe-Mass-Extin ction distribution ob- 
tained with the method of Fouesneau fc LanconI 
Moise is added to the dot positions along 
the different directions to reduce the overlap in- 
troduced by the binning method. A dispersion of 
0.15 dex in age is used for this purpose, which cor- 
responds to the minimal resolution determined in 
Sect.m A dispersion of 0.05 dex is used in mass 
and 0.1 in Ay. Colors refer to the third esti- 
mated parameter, according to the scale on the 
right-hand side. 
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Fig. 6. — Age distribution obtained based on dis- 
crete population models. The top panel shows the 
dN/d\og{A) distribution whereas the bottom one 
shows the \og{dN/dA) distribution, which are the 
two commonly used representations in the litera- 
ture. The black curves represent the distribution 
for the whole set of detected clusters, whereas the 
gray curves only include objects with estimated 
masses above lO'^'^ Mq. Error bars are the Pois- 
son noise dispersions in each bin. Solid and dashed 
lines represent power-law fits to the distributions 
of the bottom panels (for log(^/yr) between 7.0 
and 9.0). The bin width is of 0.35 dex. 
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ters come with a large range of extinction val- 
ues, as seen in many sta r forming galaxies (e.g. 
Whitmore fc Zhang|[20ol Kim et al. in prep.). At 
ages older than 10'' yr, clusters with more than 
one magnitude of extinction become rare. Be- 
tween 10'' and 10^ yr, this most likely reflects 
a real lack of highly reddened objects, as the 
data's detection limits in principle allow us to 
detect 3 x 10^ yr old clusters with masses above 
log(M) ^ 3.7 and up to 2 magnitudes of extinc- 
tion. At old ages, we do not expect many highly 
reddened clusters to be present in the disk of M 83. 
But even if they existed they would have faded be- 
low the detection limits of the sample unless they 
were very massive. 

The bottom panel is consistent with expecta- 
tions again: at a given age, the upper envelope of 
the derived extinctions decreases with decreasing 



6.1. Cluster ages 

The distribution of all the cluster ages is shown 
in Figure El The age distribution of the clus- 
ters can be approximated by a power law, with 
an index —0.93 ± 0.10 for clusters more massive 
than log(M/ Mq) = 3.5 and with ages between 
w 10'^ - 10^ yr. If we fit from 3 x 10^ to 10^ yr 
instead, we find a power law index of— 1.12±0.18. 
The slope remains similar if we restrict the sample 
to even higher m asses. These results are similar 
to those found by I C bandar et al. ( 2010[ ). although 
they used different ranges of masses and signifi- 
cantly broader bins to account for age-dating ar- 
tifacts resulting from continuous models. 

For ages between 10'' and 10^ yr, the sample's 
age distribution is roughly flat in logarithmic age 
bins (top panel). Substructures are observed in 
some age distributions with different bin sizes and 
locations. These substructures disappear when 
the age distribution is plotted with bins wider than 
0.4 dex. The distribution drops off at ages older 
than 10^ yr, as expected from the observational 



selection limits. 

The age distribution of the full sample has a 
peak at ages of 3 — 10 x 10^ yr. The excess of 
clusters found at young ages, compared to a power 
law extrapolation of the distribution at older ages, 
consists of low mass objects (log(M/ Mq) < 3.5). 
Objects with masses this small typically have 



faded below our detections limits when older than 
10^ yr. Imposing a mass cut at log(M/ Mq) = 3.5 
largely removes this peak. 

Systematic errors are likely to affect the strength, 
width and position of the peak at young ages. In- 
deed, the derived distribution falls to nearly zero 
below ages of about 3 Myr, despite the fact that 
M83 is still in the process of forming stars and 
clusters. It is likely that some of the clusters with 
estimated ages near 5 Myr are in fact younger. 
Reasons for possible systematics in this age range 
include: (i) the sensitivity of the age dating to 
the assumed escape fraction of ionizing photons 
from the region in which the fluxes are measured 
(as discussed in SectH]); (ii) some level of incon- 
sistency between the stellar evolution models and 
the clusters in M83; (iii) the fact that our cur- 
rent model collection contains only models with 
ages that are integer multiples of 1 Myr; (iv) the 
lower-mass limit of the current collection; and (v) 
our restriction to Ay < 3 in the analysis of cluster 
colors. A full usage of the posterior probability 
distributions of each individual cluster might re- 
duce some of these systematics (see Sect. 17. 4|) . 

The spatial distributions of young and old clus- 
ters highlight the effect of spiral arms (Fig. [71). As 
illustrated by Fig. |8l a spatial plot of the subsam- 
ple of clusters with assigned ages younger than 
10^ yr clearly shows the spiral pattern and high- 
lights the recent star formation in the nuclear re- 
gion. The spiral pattern, as traced by the clusters, 
starts fading away for clusters older than approx- 
imately 300 Myr, and the (only) 64 clusters with 
ages older than 10^ yr are consistent with a uni- 
form distribution throughout the disk. As selec- 
tion effects may differ in crowded and reddened 
areas of the field, these uncorrected spatial distri- 
butions should be taken with caution. 

In summary, the age distribution of clusters 
younger than 10^ yr and with masses log(M/ Mq) 
> 3.5 is approximately a power law with an index 
of —1.0 ±0.2. While this is similar to previous re- 
sults, the distribution is now defined with better 
age resolution (0.4 vs. 0.7 dex), and it is based on 
more appropriate modelling than previous deter- 
minations when masses below log(M/ Mq) of 4.5 
are included. 
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Fig. 8. — This figure represents the spatial-age distribution of the cluster sample. On each panel, gray dots 
show the whole sample while black dots are objects with ages indicated in the top right corner. 
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6.2. Cluster masses 
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Fig. 7. — From top to bottom, recovered ages, 
masses, and extinctions vs. location in the galaxy. 
Colors are coded according to the respective scale 
and parameter given on the top of each figure. 



The marginal distribution of ah the clus- 
ter masses of the sample (i.e. the distribution 
summed over all ages and extinctions) is shown 
in Fig.ini It can be approximated by a power-law 
with index —2.15 ± 0.14 for masses higher than 
log(M/ M0) of 3.5. The slope remains within 
the uncertainties (—2.09 ± 0.13) when we re- 
strict the sample to clusters with ages between 
10^ — 10^ yr. This confirms results from previous 
analyses, which al so found approximatively power- 
law distributions (* Anders et al ■l2007tl iDowell et al 



2008; Fall ct al. 2009, 20051: IChandar et al.ll20l'or 
Bastian fc de MinkI |2009|) . Note that the power- 
law of the prior mass distribution in our analysis 
is —1: the inversion is able to move away from the 
prior if necessary. 

Below log(M/ M0) — 3.5 the derived mass 
distribution falls off, and ends with an appar- 
ent peak near log(M/ M©) = 3. The peak at 
these low masses is unreliable, as it corresponds 
to the lower mass limit of the collection of syn- 
thetic clusters used in the analysis: some of these 
objects may in fact be even less massive than 
10'^ Mq. The relative lack of clusters found be- 
tween log(M/ M0) = 3.2 and 3.5 on the other 
hand is not related to any known artifact of the 
analysis (see however Sect. 17. 4|) . It may well result 
from the incompleteness of the cluster sample at 
low fluxes: based on the Monte-Carlo collection, 
we know that log(A// Mq) = 3.7 is the threshold 
mass under which more than half of the clusters 
have fluxes below the "completeness" limit (de- 
fined in Sect.H]) in at least one photometric pass- 
band. 

Also, one should recall that (i) being faint, the 
low mass objects have larger observational errors 
than massive ones, (ii) having low masses while 
remaining above the detection limits these sources 
must be young, and therefore we are in the regime 
most sensitive to the modelling of H a. 

In summary, the use of an approach that explic- 
itly accounts for the stochastic nature of the ob- 
served clusters confirms a power-law behaviour for 
the mass distribution of the M 83 clusters, with an 
index of approximately — 2.1 ± 0.2, over the range 
of masses least affected by incompleteness, i.e. be- 
tween log(M/ Mq) = 3.5 and log(M/ Mq) = 5. 
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6.3. Nucleus versus disk 
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Fig. 9. — Mass distribution obtained with the 
stochastic analysis. This is the mass counterpart 
of the age distributions given in Fig.[6l The black 
curves represent the distribution for the whole set 
of identified objects whereas the gray curves only 
includes objects with estimated ages higher than 
10 Myr. Error bars are the Poisson noise disper- 
sions in each bin. Power-law fits are overlaid on 
top of the distributions (fitted masses are above 
lO^-s Mo). 



The ^ 200 star clusters of our sample that 
are part of the galaxy nucleus have derived ages 
younger than 10 Myr with our analysis 4 Myr 
on average). This statement is not to say the nu- 
cleus contains no older objects, since most older 
clusters would remain undetected against the high 
background of the young sources. The higher av- 
erage mass of the nuclear clusters ^q) is 
also compatible with the selection effects expected 
in an area crowded with bright objects. 

Our main conclusions about the age distribu- 
tion and the mass distribution of the M 83 sam- 
ple remain valid for the disk when the nuclear 
subsample is removed. The age distribution is 
changed only below 10 Myr, hence the fits shown 
in Fig. [6] are essentially unaffected. The field clus- 
ters' mass distribution is only slightly depleted at 
high masses compared to the mass distribution of 
the whole sample, and this steepens the index of 
the fitted power law by only a few percent. 



7. Discussion 

7.1. The importance of the 
F657N 



Ha filter, 



Despite the importance of Ha in the identifica- 
tion and age-dating of young clusters, this filter is 
still not used extensively. This is partly because 
extra difficulties come with these measurements 
and their interpretation. For example, the spatial 
extent of HIT regions rarely matches that of the 
underlying optical continuum light leading to dif- 
ficult aperture corrections, the escape fraction of 
ionizing photons is poorly known, the hot ionizing 
stars are a rare subpopulation and their number 
in any given clu ster can deviate stro ngly from the 
expected mean (jCervino et al.ll2003l ). Considering 
these difficulties, we tested the results of the anal- 
ysis of the M83 clusters with UBVI data alone, 
i.e. discarding the flux measurement in the Ha 
filter (F657N). 

Figure [10] shows how the ages and masses es- 
timated with Ha included (y-axis) compare with 
those obtained from UBVI alone (x-axis). In this 
figure, the colors represent changes in the estimate 
of the extinction parameter Ay between the two 
studies. The age-extinction degeneracy is appar- 
ent: younger ages come with naturally bluer col- 
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Fig. 10. — Comparing age and mass esti- 
mates (from the stochastic approach) based on 
UBVI+Ha (y-axes) and UBVI only (x-axes). 
Each dot represents the estimates of a star clus- 
ter color-coded according to the change in the 
extinction estimate (resulting from the inclusion 
of Ha). Age estimates are compared on the top 
panel, while mass estimates are compared on the 
bottom panel. Projected age and mass distribu- 
tions are represented in Fig. 1111 
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Fig. 11. — Comparing age and mass distribu- 
tions (from the stochastic approach) based on 
UBVI+Ha (blue) and UBVI only (black). Age 
estimates are compared on the top panel, while 
mass estimates are compared on the bottom one. 
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ors on average, which can be compensated with a 
higher Ay. The blue dots show that without the 
narrow-band filter information many clusters with 
Ha emission are given intermediate ages (40%). 
The masses of the corresponding clusters are sys- 
tematically lowered when the Ha filter is included: 
the extra extinction required to match the colors 
at the now younger ages does not quite compen- 
sate for fading (younger clusters are intrinsically 
brighter for a given mass). The red dots also illus- 
trate the importance of R-band information. In- 
deed when there is no Ha emission, the F657N fil- 
ter becomes essentially a narrow R-band filter and 
acts as an extra constraint on the ages of clusters 
too old to contain ionizing stars. 

Figure [TT] presents the projected age and mass 
distributions. Clearly, the age distribution de- 
pends strongly on whether or not Ha narrow-band 
(F657N) data is available. However the effect on 
the mass distribution is limited: in particular, we 
conserve a similar power-law index. This suggests 
that the mass distribution of the sample based on 
our probabilistic approach is determined robustly 
down to about log(M/ Mq) — 3.5. 

The age distributions with and without F657N 
differ, but which one is closer to reality? How sen- 
sitive are the young ages to details in the modelling 
or errors in the measurements? A color-color dia- 
gram involving F657N helps answering this ques- 
tion (Fig. [H]). The locus of the cloud of data 
points is extremely well represented by the collec- 
tion of stochastic models with the extra degree of 
freedom provided by the reddening vector (this is 
reflected in the very good best-^^ values obtained 
with the stochastic models). As opposed to the 
continuous models, the stochastic ones reproduce 
the "corner" of the observed distribution at (V- 
I)~0 and (Ha -V)~ 0. 

Errors would have to be unrealistically large 
and systematic to avoid an interpretation of 
Fig. [12] with a large fraction of objects younger 
than 10 Myr. Changing the escape fraction of 
ionizing photons in the models would modify only 
the extent and color of the "plume" of young mod- 
els in the figure. Changing the slope of the upper 
stellar IMF, the exact ages for which models are 
computed, or the prior age and mass distributions 
of the model collection, would change how mod- 
els are distributed along this plume and near the 
blue "hook" in the distribution. The slope of the 
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Fig. 12. — Models and M 83 cluster data in a color- 
color plane that highlights the role of F657N. The 
model distribution is shown as a 2D histogram in 
the top panel, with the line of continuous models 
overlaid. The models are shown as a colored con- 
tour maps in the bottom panel, together with the 
M 83 data. 
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extinction vector determines where exactly the 
dereddening hne of a cluster reaches the cloud of 
models. But none of these changes would affect 
the fraction of young clusters greatly. 

To summarize, the conclusion that about 55 % 
of the clusters in the sample are less than 10^ yr 
old is robust. On the other hand, individual ages 
below 10^ yr are sensitive to model details. 

7.2. Comparison with the results based on 
continuous models 
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Fig. 13. — Age-mass distribution obtained from a 
X^-fit using continuous models. Colors other than 
cyan and green highlight clusters for which the 
ages have changed most strongly between the con- 
tinuous and the stochastic analysis (actually the 
scale is set by the corresponding change in Ay). 

For a direct comparison, we have repeated the 
analysis of the M 83 clusters with the traditional 
approach, based on the mean fluxes predicted 
by continuous models. The population synthe- 
sis assumptions are those of Pegase as described 
above, and the calculation used to measure the 
quality of a fit is de fined using fluxes as in the 
stochastic context (see Fouesneau &: Lan(;;onll2010l 
who present such a comparison for synthetic clus- 
ters). 

The 2D age-mass distribution resulting from 
the "continuous" analysis is shown in Fig.lT^ 
When compared to the results of the stochastic 
analysis (top panel of Fig.[S|) it sho ws a less homo- 
geneo us distribution in age. Foues neau &: Lanqon 
( 20101 ) explained why this is expected. 



The "continuous" marginal distributions of 
cluster masses and cluster ages are compared to 
the stochastic ones in Fig.[T31 As already stated 
in the previous section, we consider the flatness 
of the "stochastic" age distribution between 10^ 
and 10^ yr (in logarithmic bins) a robust result, 
while the dip in the "continuous" age distribution 
around 5 lO'^ yr and the excessive accumulation of 
very young clusters are artifacts of that method. 

To first order, the mass distributions in the 
bottom right panel of Fig. [14] are alike. Above 
log(M/ M0) = 4, they are within 2 ct of each other 
(ct representing Poisson standard deviations in the 
bins). 

A closer look shows that the analysis with con- 
tinuous models leads to a subset of very young 
and low mass clusters (log(M/ Mq) < 3) that are 
not found with the "stochastic" analysis. At first 
glance, one may think this difference simply re- 
flects a bias of the "stochastic" analysis, due to the 
lower mass limit of the model collection (10'^ Mq). 
In fact, it is due at least in part to a known issue 
of the a nalysis of finite clusters wi th continuous 
models ( Fouesneau fc Langon 2010l their Fig. 9). 
A purposely designed test was run with the WFC3 
filters of the M83 observations. The analysis of 
synthetic clusters, all more massive than 10^ Mq 
and brighter than the magnitude limits of the 
M83 observations, shows that '^35% of the clus- 
ters with actual log(M/ M©) between 3.2 and 3.7 
are assigned masses below log(M/ Mq) = 3 when 
analysed with continuous models. In the upper 
right panel of Fig.[T3] we see that the difference 
between the distributions derived with stochastic 
and continuous models for 3.2 < log(M/ Mq) < 
3.7 represents just about ~ 35 % of all the clusters 
in that mass range. 

More details on the differences between the 
stochastic and the continuous analysis are given 
in Appendix A. In particular, the one-by-one com- 
parison of the "continuous" and the "stochastic" 
ages and masses are shown. Individual ages differ 
by more than 0.3 dex for 28 % of the sample. 

A comparison between the "continuous" results 
found in the present study (i.e. based on Pegase 
m odels analog to t he sto chastic models) and those 
of IChandar et al.l (|2010l ) is provided in Appendix 
B. Using different population synthesis models 
changes the ages of 25 % of the clusters by more 
than 0.3 dex, but in a way that does not modify 
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Fig. 14. — Marginal distributions of cluster ages and masses as obtained with continuous models (solid 
and dashed histograms) and with stochastic models (red histogram levels). In each case, the histograms are 
shown for two sets of bin centers, offset from each other by half a bin size. To guide the eye, error bars 
representing Poisson standard deviations are shown in the top panels. 



global results. 

7.3. Comparison w^ith morphological or 
spectroscopic ages 

Age estimates based on methods other than the 
analysis of the broad and narrow-band Spectral 
Energy Distributions (SED) are available for small 
subsets of the star clusters of our sample. Unfor- 
tunately these subsets do not extend to the lowest 
cluster masses. 



Whitmore et al 



()201lD have devised a purely 
morphological classification scheme for star clus- 
ters based on images in several photometric pass- 
bands, including Ha. The categories are chosen 
such as to represent subsequent phases of the mor- 



phological evolution of a typical cluster, and the 
classification sequence is therefore expected to cor- 
relate with age. Their Fig. 3 shows that the mor- 
phological category indeed relates to the SED ages 
estimated using traditional models with a contin- 
uously populated IMF. 

^_Fig JT5l focus e s on the 64 star clusters of 
Iwhitmore et al. ( 2011 ) for which the morphologi- 



cal classification is considered most robust. Eight 
of the 64 clusters are located in the nucleus, 56 
in the field. The bottom panel compares our age 
estimates with the pho tometric ages obtained by 



Whitmore et al.l (|2011l ) based on continuous mod- 



els. Both of these age estimates use UBVI-|-Ha 
photometry. The two sets of photometric ages 
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Fig. 15. — Ages derived from the approach 
detailed in this paper [a;-axis] c ompared to the 



Ha m orphology categories from iWhitmore et al 



( 201l[ ) (top) and to ages from traditional models 
(bottom), for the clusters in common. Boxes on 
the top panel refer to statistics on the morpholog- 
ical categories presented in the cited paper. 



agree well, except for a scarcely populated hori- 
zontal branch of "catastr ophic differences" around 
\og{A/yr) = 6.9 from IWhitmore et all (|201ll ) 
(only 10% of the sample) and an associated gap in 
the one-to-one relation around \og{A/yr) = 7.2, 
a behaviour r ather ty pica l for this sort of plot 
( Fouesneau fc LanconI l201(i ). The deviant clus- 



ters are in the morphological c a tegor ies 5a and 
5b according to Whitmore et al. ( 201l[ ). i.e. they 
are unlikely to be as young as the photometric 
analysis with continuous models suggests. 

The top panel of Fig. [15] shows the relationship 
between photometric age and morphological cate- 
gory obtained with the new photometric age esti- 
mates, i.e. those based on stochastic models. The 
clusters that are outliers in the upper panel were 
also outliers in the ver sion of this figure given by 



Whitmore et al.l (|2011i ): they formed a vertical ex- 



tension of young photometric ages into older mor- 
phological categories. But for one object, this ver- 
tical extension has now disappeared. As a result, 
the average photometric ages of categories 5a and 
5b shift to older values by about 0.1 dex. This re- 
mains within the error bars associated with the 
dispersion within the categories. All in all, the 
figure remains similar to its original version. 

Spectroscopic age estimates are available for 
13 luminous c l usters in the nuclear region from 
Wofford et al.l ( 2011 ). These were obtained us- 



ing STIS FUV spectra covering the range from 
1200 - 1700 A (e.g., including the strong N V 
and C I V lines as wel l as sev eral other weaker 
lines). IWofford et al.l (|201ll ) fit the observa- 
tions using both semi-empirical models, based 
on a library of Galactic O and B stars observed 
with lUE, and th e oretic al Starburst 99 models 
(|Leitherer fc ChenI l2009l continuously sampled 
IMF). They find ages between 2 and 5 Myr for 
nine clusters and ages between 10 and 20 Myr for 
four. The photometric ages quoted in that ar- 
ticle, based on the same five photometric pass- 
bands as in our pa per and on the method of 
Chandar etaP (|2010l ). are in good agreement with 



the spectroscopic ages except for one of the four 
obj ects with spectro scopic ages above 10 Myr (see 
I Wofford et al.ll201l[ Table 3). Our analysis with 
stochastic models returns ages between 3 and 
9 Myr for all the 13 clusters. These ages are within 
0.2 dex of the spectros copic results except f or the 
three oldest objects of I Wofford et al.l (|201l[ ). The 
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differences for these three objects are due to differ- 
ing assumptions in the population synthesis codes 
rather than to stochastic modelhng. Indeed, when 
we use models with a continuously sampled IMF 
we obtain ages within 0.1 dex of those derived in 
the stochastic context, as long as both are based 
on the same synthesis parameters (e.g. those de- 
fined in Sect. [3]). This comparison is meaningful 
because the three clusters are massive (masses in 
the range 5 x iC — 1 x 10^ Mq depending on the 
method used). All in all, the agreement between 
stochastic and non-stochastic ages is very satisfac- 
tory considering the diversity of approaches. 

7.4. Towards a more complete usage of 
posterior probability distributions 

The main aim of this paper is the comparison 
of results of traditional anal ysis methods with re- 
sults b ased on the method of lFouesneau fc Lanqon 
( 2010l ). for a sample of real clusters in M 83. Fol- 
lowing these authors, we have chosen the peak of 
the posterior probability distribution to assign an 
age, a mass and an extinction value to each cluster. 
This method can be extended to use the full poste- 
rior probability distribution for each cluster rather 
than the peak approximation. In this section, we 
explore this approach by comparing the results 
of this extended method with those of the peak 
method for our cluster sample in M83. A com- 
plete theoretical study of the extended method in 
various observational regimes, will deserve a sepa- 
rate paper. 

The 3-dimensional probability distributions of 
age, mass and extinction obtained for the 1242 in- 
dividual clusters of our sample are quite complex 
and display a variety of patterns. We have found 
about a dozen typical behaviours, that depend on 
the region of color-luminosity space the clusters 
are in, as well as on observational error bars. It 
is difficult to find a satisfactory way of translat- 
ing this complexity into simple error bars on the 

estimates provided by the "peak" metho d. As 

already described in lFouesneau fc LangonI ( 2010l ). 
the main trends are set by the age-mass and age- 
extinction degeneracies, which define directions 
along which probability peaks tend to be elon- 
gated. Within this gross picture other patterns are 
seen. For instance, it is relatively common to find 
bimodal distributions with one peak at young ages 
(< lO^r) and one at intermediate ages(~ 10^ yr). 



The relative strengths of these peaks depend on 
the error bars on the Ha measurement. This par- 
ticular pattern can be understood with the discus- 
sion we provide in Appendix A concerning the left 
panel of Fig.Ull 

A new estimate of the age and mass distribu- 
tions of the cluster sample as a whole was obtained 
by summing the age-mass probability distribution 
of individual clusters (marginalized over extinc- 
tion). With this approach, uncertainties on the 
estimated parameters and any covariance between 
them are automatically accounted for. The sensi- 
tivity to binning choices is reduced. The resulting 
map is shown in Fig. 1161 For convenience, the fig- 
ure also recalls the results of the "peak" method 
(Fig.[5|) as a directly comparable density map. The 
global structure with both approaches is similar. 
For instance, both methods obtain an over-density 
of young, low-mass objects. Overall, by account- 
ing for the entire posterior probability distribu- 
tions instead of only the location of their maxima 
we smear out local variations in the derived age- 
mass map. 

The one-dimensional age and mass distribu- 
tions obtained with the extended approach are 
presented in Fig. ll7[ and can be compared to those 
shown in Figs. [6] and [9l The masses of the cluster 
sample are found to follow a power-law with an 
index of —2.1 for log(M/ M0) > 3.5, a result es- 
sentially identical to the one we obtained with the 
peak method. The restriction to ages older than 
10 Myr does not change the index value signifi- 
cantly. The dip seen near log(Af /M©) = 3.3 with 
the peak method has been smoothed into a more 
robust plateau. The age distribution in Fig.fTTl 
has a somewhat steeper power-law (a = —1.23) 
than found previously (—1.08). Its restriction to 
masses greater than log(Af/M0) = 3.5 reduces 
the absolute value of the age power-law index by 
some 10%. The slopes of the age distributions, 
as derived with either the extended and the peak 
method, remain within the formal 1-a confidence 
interval of each other for this sample. 

In summary, the age-mass distributions found 
for our sample of M 83 clusters are similar whether 
we use the extended method or the peak method. 
Spurious small scale features in the result of the 
peak method tend to disappear with the extended 
method. The extended method deserves further 
investigation in the future. 
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Fig. 17. — Age and mass distributions obtained with the approach explored in Sect. 17. 41 Left: age distribu- 
tion, for comparison with Fig. [51 Right: mass distribution, for comparison with Fig. [HI 



21 



8. Conclusions 
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Fig. 16. — Age-mass distributions obtained for 
the M 83 sample from the probabilistic analysis. 
Top: distribution obtained by summing the indi- 
vidual posterior probability maps of the clusters 
(i.e. the extended method). Bottom: density 
map obtained from the most probable age and 
mass values shown in the top panel of Fig. 5 (i.e 
the peak method, for which each point is replaced 
by a patch to produce a density representation). 



Most star clusters in nearby galaxies have 
masses below 10^ M©, and hence their inte- 
grated fluxes and colors can be strongly af- 
fected by stochastic fluctuations in the num- 
ber of massive stars (strictly post main-sequence 
stars). iFouesneau fc LanyonI (1201 Q ) previously in- 
troduced new stochastically sampled stellar pop- 
ulation model predictions, and developed a prob- 
abilistic method for estimating the ages, extinc- 
tions, and masses of star clusters from these mod- 
els. The models were tested using synthetic clus- 
ters. For the purpose of the present paper we 
enhanced the stochastic models by including pre- 
dictions for narrow-band filters and extending the 
Monte-Carlo collection of synthetic clusters to in- 
clude higher cluster masses. We compared the pre- 
dictions from the stochastically sampled models 
with observations of a sample of real star clusters. 
Our observations consist of integrated UBVIHa 
fluxes of star clusters in the nearby spiral galaxy 
M83, which was observed with the Wide Field 
Camera 3 (WFC3) on the Hubble Space Telescope. 
This is the same cat alog of clusters analysed by 
Chandar et all (|2010t ) using predictions from more 
traditional models which assume a "continuously 
sampled" stellar IMF. 

The locus of the collection of stochastic mod- 
els in color space (e.g.. Figure [5]) shows excellent 
agreement with that of the collection of cluster ob- 
servations. Clusters with broad band colors either 
bluer or redder than those of the traditional mod- 
els find a natural match in the stochastic model 
collection, in which some synthetic clusters have 
underpopulated or overpopulated red post-main 
sequence branches. Similarly, the locus of stochas- 
tic models in color-color planes that include Ha 
fluxes (e.g.. Figure [T2|), and hence take into ac- 
count the random character of the number of ion- 
izing stars in a cluster, provide a natural match 
with the observations. 

A key result of the study is the importance of 
including a narrow band filter measurement in the 
analysis, which to date has been used in relatively 
few works. A comparison of age and mass esti- 
mates with and without the Ha filter shows that 
30% of the sources have age estimates that change 
by more than 0.3 dex in age in our analysis (factor 
of 2). The errors are systematic: without the nar- 
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row band measurement, many young clusters are 
not recognized as such. 

The use of stochastic models has allowed us 
to derive the age-mass distribution of the star 
clusters in the M 83 sample with a better reso- 
lution in time than was possible before, and we 
extend both age and mass distributions to lower 
masses (down to log(M/ Mq) ^ 3.3). Above 
log(M/ M0) = 3.7, the mass limit at which com- 
pleteness issues in the sample become severe, the 
mass distribution derived from the analysis with 
stochastic models closely resembles the distribu- 
tion obtained with traditional models in previous 
work. More specifically, the overall slope of the 
mass distributio n remains within the range of val- 
ues indicated bv lChandar et al. I (l2010l) . I.e. power- 
law of index near —2. The age distribution is now 
free of the dip between 10^ and 5 x 10^ years that 
was a known artifact of the analysis of clusters 
subject to stochastic fluctuations with continuous 
models. Overall, we confirm that it declines more- 
or-less continuously starting at yo ung ages as al 



ready suggested by the results of IChandar et al 



( 2010l ). Between 10^ and 10^ years, the data are 



consistent with a flat distributions in logarithmic 
age units, which corresponds to dN/dA oc A~^. 
The sample contains only a handful of clusters 
older than 1 Gyr because these are too faint for 
easy detection. 

In this paper we also begin the process of ex- 
tending the analysis from a unique "most proba- 
ble" value (corresponding to the peak in the poste- 
rior distribution, the peak approach) to a full prob- 
abilistic description of the correlated age, mass, 
and extinction properties of each individual clus- 
ters (i.e. the extended method). The extended 
method is expected to be more robust than the 
peak method. It does not loose any of the infor- 
mation carried by individual posterior probability 
distributions about uncertainties, about correla- 
tions between parameters, or about multiple (al- 
most) equally probable solutions. The age-mass 
distributions obtained with the extended method 
for our M 83 cluster sample remain similar to those 
obtained previously. Spurious substructure seen 
in peak method distributions when plotted with 
narrow bins disappears. We plan to test the ex- 
tended method more systematically in various ob- 
servational regimes in the future. 

It may seem surprising at first glance that low 



resolution age and mass distributions of cluster 
samples are not modified more radically by the 
move from a traditional analysis to a method 
based on explicitly stochastic models. Changing 
a power law slope radically, with data that ex- 
tends over about two dex in age and mass, re- 
quires strong and systematic changes, but the 
changes induced by the switch from traditional 
model s to stochastic mode l s are not of this na- 
ture. Fouesneau fc Lanpon ( 2010l) had predicted 
that overall mass distributions would be rather 
insensitive to this switch, and we have conflrmed 
this with the M 83 sample. Their figures also indi- 
cated that changes in the age assignments due to 
the switch to stochastic models would mostly af- 
fect subsets of objects in particular ranges of age 
and mass (i.e. the high resolution features of the 
age distribution), instead of changing the global 
balance between young, intermediate and old ob- 
jects. This is also what we found for M 83 sample 
in the present paper: the —1 (±0.2) slope of the 
best fit power law to the age distribution is a ro- 
bust result for the sample. 

As we push to lower masses and larger sam- 
ples of fainter clusters, the improved accuracy and 
time resolution achievable with the new stochas- 
tic methods allows us to address new questions, 
such as local variations among cluster populations 
within individual galaxies. For instance, we find 
that the subsample of young clusters in our data 
set clearly displays the spiral structure of M83, 
while this structure progressively disappears in 
samples older than a few hundred millions of years. 
The cluster age-mass distribution in the field we 
have studied is thus (not surprisingly) partly de- 
termined by the location of this field with respect 
to the spiral arms. This level of details can now 
be taken into account when comparing galaxies or 
areas within galaxies. 

This paper is based on observations taken with 
the NASA/ESA Hubble Space Telescope obtained 
at the Space Telescope Science Institute, which is 
operated by AURA, Inc., under NASA contract 
NAS5-26555. It uses Early Release Science obser- 
vations made by WFC3 Science Oversight Com- 
mittee. The analysis of the results made extensive 
usage of the Topcat software, available under Gen- 
eral Public Licence from 



http : //www ■ starlink . ac . uk/topcat/[ 
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A. Appendix A: More details on the comparison between the continuous and the stochastic 
analysis 

This appendix provides some additional details about the comparison described in Sect. l7.2l 
In Fig. llSi we provide the direct chister-by-cluster comparison between the ages and masses derived on 
one hand from an analysis with stochastic models, on the other hand from the best-x^ fit to continuous 
models (both based on the population synthesis code Pegase, and on UBVI and Ha photometry). It is 
essential to note that in this section (as in most of the paper) the parameters we assign in the stochastic 
context are not those of the single Monte-Carlo model that provides the best fit, but those that maximize 
the posterior probability distribution of the cluster observations. 72 % of the clusters have ages within 0.3 dex 
of each other, 28 % have ages more than 0.3 dex apart. For the clusters with similar ages in both approaches, 
the residuals have a bell-shaped distribution with a FWHM of 0.3 dex. Most of the clusters with highly 
method-dependent ages lie on a vertical sequence in the second panel of Fig.[TH] (21 % of all clusters, see 
the box on the figure). They are assigned young ages with continuous models, while the analysis based 
on stochastic models leads to a range of older ages. The young ages of the continuous analysis come with 
high values of the extinction parameter (because young cluster models are intrinsically blue) and low masses 
(because young cluster models are intrinsically bright), while the older ages of the stochastic analysis are 
associated with low extinction parameters, and masses that are typically 0.8 dex higher. This is seen in 
the first panel of the figure as a one-sided broadening of the one-to-one relation towards higher "stochastic" 
masses (dark blue symbols). Because this happens at all masses, the consequence for the mass distribution is 
a slight offset in log(Af/ Mq) with hardly any change in the slope of a fitted power law. The consequence for 
the age distribution is a change in the strength of the peak of young clusters (< 10^ yr), and a redistribution 
of these clusters over other ages in a way that again has little effect on the slope of a fitted power law. 
Note that for comparison with other work in the literature one also has to keep in mind differences in the 
population synthesis assumptions of various authors (see e.g. Appendix B). 

The remainder of this appendix is provided for readers particularly interested in the subtleties of the 
behaviour of our analysis with stochastic models. It gives further insight into the roles played by the prior 
distributions of cluster ages and masses, and by the observational error bars. 

Color-color and color-magnitude plots can be used to understand the origin of the vertical sequence in 
the age diagram just described. For brevity we will call these particular clusters the "deviant clusters" in 
this paragraph. We looked at many color combinations and confirmed that the deviant clusters are located 
in a region of broad band color space where young reddened models and old models overlap. The projection 
we found most useful for this discussion is (Ha —V) vs (V— I). The deviant clusters are highlighted in this 
diagram in the first panel of Fig. [HI which must be compared to Fig. [121 About one half of the deviant 
clusters (shown in blue and cyan) lie along the line where most of the old and intermediate age synthetic 
clusters congregate, near the upper envelope of the total cluster sample in this diagram. The fact that they are 
assigned intermediate and old ages in the stochastic analysis is an immediate consequence of the high density 
of intermediate and old model clusters in that area. The fact that some are assigned very young ages in the 
continuous analysis appears as a consequence of the "single best fit" approach of that continuous analysis: 
by chance the model closest to them happens to be a young one, but a small shift in the observed colors 
(well within the photometric error bars) could change the derived age significantly since young reddened 
models and older dust-free ones are found side by side. The second half of the deviant clusters (shown 
in yellow and red) are located in a region where only young reddened model clusters exist. These objects 
however all have errors larger than 0.25 magnitudes on the F657N measurement (and the errors increase 
for deviant clusters located further away from the upper envelope of the data points). For these clusters 
we see the effects of the observational errors clearly. The multi-dimensional 2 a error boxes around each of 
them reach well into the above-mentioned region of color-color space where the density of intermediate age 
and old clusters is high. Models outside such an error box contribute negligibly to the posterior probability 
distribution, but all models inside may contribute significantly (Bayes' theorem). Because the error boxes 
contain overwhelmingly more intermediate and old models than young reddened models, our analysis favours 
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Fig. 18. — Comparison of the ages and masses obtained when analysing the M83 cluster sample with 
"continuous" models on one hand and with "stochastic" models on the other. The dashed line highlights 
the one-to-one relation, and the dotted box the vertical sequence mentioned in the text. 
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Fig. 19. — Color-color plots similar to Fig.[T21 but with specific objects highlighted. In the first panel, 
the colored symbols identify the clusters that populate the main vertical structure of the second panel of 
Fig. [18] (i.e. the "deviant clusters"). The color code shows the observational uncertainty on the clusters' Ha 
measurement (in magnitudes). In the second panel, the colored symbols highlight low mass clusters whose 
ages are similar in the stochastic and in the continuous contexts, but whose masses differ. Here the color 
code shows the associated Ay. Both subsets are discussed in the text. 
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the older solution. For clusters at the same location in color-color space but with much smaller error bars 
on Ha , young ages are returned. Note that the number of clusters concerned by this discussion is too small 
to change any of the main results discussed in this paper. 

A second more subtle effect contributes to the one-sided broadening of the mass vs. mass plot, this 
time only in the regime of small masses (log(Af/ Mq) < 3.5 in the continuous analysis - red symbols in 
the first panel of Fig. lTS)) . These objects have identical ages with the continuous and with the stochastic 
analysis (frequently within 0.1 dex), but the "stochastic" masses are larger than the "continuous" masses 
by 0.3 — 0.6 dex. Their location in the most useful color-color diagram is highlighted in the second panel of 
Fig. [121 We see in Fig. [T^] that the vast majority of the young stochastic cluster models lie not along the line 
of models with a continuously sampled IMF, but rather on a curved "plume" at bluer colors. This would 
be the case for any prior mass distributions that has more low mass clusters than high mass clusters. As a 
consequence, the analysis of the posterior probability distribution favors the models in this bluer "plume" 
and assigns a correspondingly larger extinction. The age changes little when compared to the analysis with 
a continuously sampled IMF, but the larger extinction must be compensated for with a somewhat larger 
mass. This effect concerns about 20 % of the clusters. It has a noticeable effect on the mass distribution 
derived from stochastic models only below log(Af/ M0) — 3.6. 



B. Appendix B : comparison with the ages and masses of Chandar et al. 2010 

Chandar et al. 1 (I2OIOI) use continuous population synthesis models. Figure 10 of their article shows the 



age-mass distributions they obtain when assuming solar and twice solar metallicity. They are qualitatively 
very similar to the distribution shown here in Fig. 1131 

The d ifferences are due in pa rt to the underlying population synthesis assumptions, in part to the analysis 
method. IChandar et al. ( 2010[ ) use models from Bruzual & Chariot (2009) with a Chabrier IMF, and they 
account for hydrogen line emission with an escape fraction of 0.4 (the value that, in their analysis, provides 
best fit qualities on average). The fit for age and extinction is done using magnitudes, and the mass is then 
determined from the V band absolute magnitude. As a consequence of using magnitudes in their case, fluxes 
in our case, the different photometric bands are given different weights in both studies, and the reduced 
values cannot be compared directly (except very grossly: in both cases a value above a few indicates a 
poor fit). 

Figure[20l directlv compares ages retrieved bv IChandar et aL ( 2010[ ) and ages we obtain in this present 
paper. Ages are similar for 85 % of the clusters (and for 95 % of the clusters with small observational errors). 
These objects are on a curved line that deviates slightly from the one-to-one diagonal, so that strictly 
speaking only 75 % of all the clusters have ages that agree to within 0.3 dex with both methods. Differences 
one could call "catastrophic" occur for 15 % of the clusters, and are not restricted to observations with large 
error bars. These differences depend rather sensitively on the way nebular emission is accounted for in the 
models (the vertical branch of "catastrophic errors" disappears and the horizontal one gains clusters when 
we remove the nebular emission from the continuous models used here). They also depend on the range over 
which the extinction parameter is allowed to vary. 
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Fig. 20. — Ages derived from with continuous models by IChandar et all (|2ni(t [ X-axis] and in this paper 
[y-axis]. Large black dots highlight the 114 clusters for which the error in all filters are smaller than 0.13 
mag. 
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